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Motivated from recent experiments in k-(BEDT-TTF)2Cu2(CN)3, dielectric and magnetic 
properties in a dimer-Mott insulator are studied. We derive the effective Hamiltonian where 
the number of electron per dimer is taken to be one. Charge and magnetic structures in finite 
temperature are obtained by using the mean-field and Monte-Carlo methods. Magnetic and 
ferroelectric phases are exclusive with each other. Dielectric fluctuation is remarkable near a 
phase boundary between dimer-Mott phase and ferroelectric phase. Implications of the present 
results for k-(BEDT-TTF)2Cu2(CN)3 are discussed. 
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Novel dielectric and magneto-dielectric phenomena are 
one of the recent central issues in solid state physics. Be- 
yond a conventional picture based on classical dipole mo- 
ments, ferroelectricity on which electronic contribution 
plays crucial roles is termed electronic ferroelectricity.^' ^' 
Recently discovered multiferroics are known as phenom- 
ena where ferroelectricity is driven by a spin ordering. In 
a Mott insulator with frustrated exchange interactions, 
an electric polarization is induced by the exchange stric- 
tion effects under a non-collinear spin structure. 

There is another class of electronic ferroelectricity; 
charge-order driven ferroelectricity where electric polar- 
ization is caused by an electronic charge order (CO) with- 
out inversion symmetry. This class of ferroelectricity is 
observed in transition-metal oxides and charge-transfer 
type organic salts, for exmple LuFe204,3-5) (TMTTF)2X 
(X: a monovalent cation),^-^) and a-(BEDT-TTF)2l3.^) 
A large magneto-dielectric coupling and fast polarization 
switching are expected in charge-driven ferroelectricity, 
since the electric polarization is governed by electrons. 

Dielectric anomaly recently discovered in a quasi-two 
dimensional organic salt k-(BEDT-TTF)2Cu2(CN)3 sug- 
gests a possibility of electronic ferroelectricty in this 
compound. The crystal structure consists of an alter- 
nate stacking of the BEDT-TTF donor layers and the 
Cu2(CN)3 acceptor layers. In a quasi-two-dimensional 
BEDT-TTF layer, pairs of dimerized molecules locate in 
an almost equilateral triangular lattice. When two dimer- 
ized molecules are considered as a unit, since average hole 
number per dimer is one, this material is identified as a 
Mott insulator. One noticeable property observed exper- 
imentally is a low temperature spin state; no evidences of 
a long-range magnetic order down to 32mK.^°'^^) A pos- 
sibility of quantum spin-liquid states is proposed. In the 
recent experiments,^^'* temperature dependence of the di- 
electric constant has a broad maximum around 25K and 
shows relaxor-like dielectric relaxation. Some anomalies 
are also seen in the lattice expansion coefficient and spe- 
cific heat around 6K.^^'^^) These data promote us to re- 
examine electronic structure in dimer Mott (DM) insu- 
lators. 

In this Letter, motivated from the recent experimental 
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results in k-(BEDT-TTF)2Cu2(CN)3, we study dielec- 
tric and magnetic properties in a DM insulating system. 
From a Hubbard-type Hamiltonian, we derive the effec- 
tive Hamiltonian where the number of electron per dimer 
is one. By using the mean-field (MF) approximation and 
the classical Monte-Carlo (MC) simulation, we examine 
spin and charge structures in finite temperature. It is 
shown that the ferroelectric and magnetic phases are ex- 
clusive with each other. A reentrant feature of the DM 
phase enhances the dielectric fluctuation near the CO 
phase. Implications of the present results for k-(BEDT- 
TTF)2Cu2(CN)3 are discussed. 

We start from the model Hamiltonian to describe the 
electronic structure in a DM insulator. A moleculer- 
dimer is regarded as a unit and is allocated at each site of 
a two-dimensional triangular lattice. An average electron 
number per dimer is assumed to be one. The Hamiltonian 
consists of the two terms as 



(1) 



The first term is for the intra-dimer part given by 

^intra = £^ cj^s'^vs ^ ^0 X/ i^ias'^ibs + ^ 



■ifis 



(2) 



where two molecules are identified by a subscript /i(= 
a,b). We introduce the electron annihilation operator 
Ci^s for molecule fi, spin s(=t,4,) at site i, and the num- 



ber operator 



We con- 



sider a level energy e, the inter-molecule electron transfer 
io(> 0) in a dimer, the intra- molecule Coulomb interac- 
tion Uq and the inter-molecule Coulomb interaction Vo in 
a dimer. In addition, the inter-dimer part in the second 
term of Eq. (1) is given by 
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where if/^ and V-^^ are the electron transfer and the 
'J ''J 
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Fig. 1. (Color online) Pseudo-spin directions in the — Q^ plane 
and electronic structures in a dimer. 

Coulomb interaction between an electron in a molecule 
/i at site i and that in a molecule fi' at site j, respectively. 
The first and second terms in "Hintcr are denoted by Jit 
and Hv, respectively. 

We briefly introduce the electronic structure in an iso- 
lated dimer. In the case where one electron occupies a 
dimer, the bonding and anti-bonding states are given 
by |/3«) = Has) + \bs))/V2 and = (|a,) - |6.))/^/2 
with energies Efj — e — Iq and Ea ~ e + to, respec- 
tively. In these bases, we introduce the electron opera- 
tor Ci^s for 7 = (a, /3) and the electron transfer integral 
between the NN molecular orbitals 7 and 7'. These 
are obtained by the unitary transformation from Ci^^ 
and t'^t . Two-electron states in a dimer are the follow- 
ing six states: the spin-triplet states {|T-|-), |Tj_), |To)} — 

{|at/3t),|a4/3i)>(|at/?4) + |ai/?t))/v^} ^it^ the energy 
Et — 2e + Vb, the spin-singlet state 15) = {\a^(3i) — 
\ai(3f))/V2 with £^5 = 2e + Uo, and the doubly- 
occupied states \D+) — Ci\a^ai) + C2\P'\-l3i) and \D-) = 
C 2\atai) - Ci\(3^Pi ) with Ed± = (4£ + C/o + Vb ± 
ViUo - Vo)'^ + 16t^)/2 and coefBcients C2/C1 = {Uq - 
Vo)/[2Ed+ - 4(e - to) - (Uo - Vb)] . The lowest eigen state 
is \D-). The effective Coulomb interaction in the lowest 
eigen state is Ueff = Ed_ — 2Ep ~ Vb + 2to in the limit 
of J7o,Vb » to. 

From the Hamiltonian in Eq. (1), we derive an effec- 
tive model in the subspace where each dimer is occupied 
by one electron. Charge structure in each dimer is repre- 
sented by the pseudo-spin (PS) operator with amplitude 

of 1/2 defined by = 5 E^' "^U'^mm'Vs with the 
Pauli matrices a. The eigen states for with the eigen 
values of 1/2 and —1/2 are the charge polarized states 
where an electron occupies the a and 6 molecules, respec- 
tively, and those for Qf with 1/2 and —1/2 correspond 
to the covalent states where an electron occupies the jS 
and a orbitals, respectively (see Fig. 1). 

Here we derive the effective Hamiltonian. We assume 
C/o,Vb >> jV^j^ : and treat the inter-dimer term 
Winter in Eq. (1) as the perturbed term. The effective 
Hamiltonian up to the orders of 0(^1) and 0(7^^) is 



given hy % = Hmtra + Hv + T~Lj. The first and second 
terms correspond to "Hintra in Eq. (2) and the second 
term of "Hintor in Eq. (3), respectively, where the doubly 
occupied states in a dimer are prohibited. By using the 
PS operator, these terms are given by 

^intra + - -2to ^ + ^ VF-.j Q^QJ, (4) 

i (ij) 

where W^- (= V^'j" + V,f - V,f ~ V^f) is the effective inter- 
dimer Coulomb interaction. This part is nothing but the 
transverse-field Ising model. 

The third term in H is the exchange-interaction term 
derived by the second-order perturbation with respect 
to Ht in Eq. (3). This is given by a sum of the terms 
classified by the two-electron states in a dimer de- 
noted by m as Hj — X^m ^j"^ where a suffix m — 
m, Ti, To, S, D+,D^}. A dominant term in Hj is H^f''' 
which has the lowest intermediate state energy. This is 
explicitly given by 

with 

+ E {Jl'Q>n + Jl^n,.Ql)- (6) 

7 = (",/3) 

We define Qf = Qf±iQf and n,„(/3) = l/2-(-f-)Qf. Ex- 
pressions of the exchange constants are given in Ref.^^^ 
Other terms in Hj are represented by similar forms with 
Eqs. (5) and (6). Spin and charge degrees are coupled 
with each other in the Hamiltonian, although the SU(2) 
symmetry is preserved only in the spin sector. This type 
of the Hamiltonian is similar to the so-called Kugel- 
Khomskii model^^' in an orbital degenerated Mott in- 
sulator, and was also proposed in study of the electronic 
state in a-NaVaOg.i^-i^) 

This Hamiltonian is analyzed by the MF approxima- 
tion and the classical MC method. In the MF calcula- 
tions, triangular lattices of 12 unit cells along (110) di- 
rection with the periodic boundary condition are used. 
We adopt the following 15 MP's, {S^'), (Q^) and (Sf^Q") 
where {p,v) — {x,y,z), in each (110) line. The multi- 
canonical MC simulations are performed in finite-size 
clusters of L sites {L < 96) with the periodic-boundary 
condition. We use 10^ MC steps to obtain histograms 
and 2 x 10^ steps for measurements. 

Finite-temperature phase diagrams obtained by the 
MF method are shown in Fig. 2. As for Wij in Eq. (4) 
and t^^^ in Eq. (3), we consider the following two cases in 
order to identify general and specific features in the re- 
sults: case I for a simple model, and case II for k-(BEDT- 
TTF)2Cu2(CN)3. In case I, we chose if/ = S^^^t {t > 0) 
and Wij — W{> 0). In case II, the three dominant f^ps 
estimated in k-(BEDT-TTF)2Cu2(CN)3 are represented 
by a parameter t as shown in the inset of Fig. 2(b), 
and Wij are considered based on the crystal structural 
data.^""^^) Similar phase diagrams are obtained in the 
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Fig. 2. (Color online) Finite-temperature phase diagrams ob- 
tained by MF method, (a) and (b) are for case I and case II, 
respectively (see text). Symbols Q:DM and Q:CO+DM represent 
the DM phase, and the coexistence phase of CO and DM states, 
respectively, and S:Para, S:120, S:3-fold I, S:3-fold II, and S:6- 
fold represent the paramagnetic phase, the 120° structure phase, 
the coUinear 3-fold spin ordered phase, the coplanar 3-fold spin 
ordered phase, and the coplanar 6-fold spin order phase, respec- 
tively. We chose W/to = 1 in (a), and Uo/to = 6, Vb/to = 4.5, 
Wi/tQ = -1, and W2/tQ = 0.07 in (b) where Wi and W2 are 
the NN inter-dimer Coulomb interactions along (100) and (110), 
respectively. The insets are schematic pictures of the inter-dimer 
transfer integrals. Vertical broken lines in (b) represent cases a), 
b), and c) (see text). 

two cases shown in Figs. 2(a) and 2(b) except for detailed 
spin structures. The CO phases reaUze in the region of 
large t. This phase is caused by Hv and the term pro- 
portional to QfQ^ in Hj. In low temperatures, magnetic 
phases are stabilized by the exchange interactions. One 
noticeable point is that when the magnetic phase ap- 
pears, the CO phase is suppressed [see around t/to — 0.6 
in Fig. 2(a) and around t/to = 0.8 in Fig. 2(b)]. As a 
result, a reentrant feature is observed in a DM phase. 
The reverse is also seen; as shown around T/tg = 0.1 
and t/to > 0.9 in Fig. 2(b), a slope of the magnetic 
transition temperature versus t curve decreases with in- 
creasing t, when the CO sets in. That is, the charge and 
magnetic phases are exclusive with each other. This is 
caused by a competition between the terms proportional 
to QfQj and QfQjSi ■ Sj in Hj; the former favors a 
ferro-type configuration of Q^, but the latter favors an 
antiferro-type one under the antiferromagnetic spin cor- 
relation {Si ■ Sj) < 0. 

Let us focus on case II in more detail. Temper- 
ature dependences of the charge susceptibility are 
presented in Fig. 3 for several values of t/tQ. We 
define the charge susceptibility along the y direc- 
tion as xc = (Ti)"^(E,eA-E»ei3)(EjeA-EjeB) 
[{QfQ]) — {Qi){Qj)] where two sublattices are iden- 
tified by A and B [see the inset of Fig. 3(b)]. We 
pay our attention to a region of t/to — 0.8 — 1. In 
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Fig. 3. (Color online) (a) Temperature dependences of the charge 
susceptibility xc along y in case a), and (b) those in cases b) 
and c) obtained by the MF method (see text). The inset of (b) 
is a schematic CO pattern in a triangular lattice. Ellipses and 
circles represent molecules and electrons, respectively. 

Fig. 2(b), there are three types of sequential phase 
changes with decreasing T; a) (Q:DM, S:Para)-j> (Q:DM, 
S:120), b) (Q:DM, S:Para)-^ (Q:CO+DM, StPara) ^ 
(Q:DM, S:120), and c) (Q:DM, S:Para)-^ (Q:CO+DM, 
S:Para) (Q:CO-|-DM, S:6-fold) where abbreviations 
are defined in the caption of Fig. 2. In the case a), al- 
though the CO phase does not realize, xc increases in 
the paramagnetic phase and suddenly reduces in the low 
temperature magnetic phase. This enhancement of xc is 
remarkable at a vicinity of the phase boundary between 
the CO and DM phases (see a curve for t/to — 0.828 in 
Fig. 3). In the cases b) and c), the CO phase realizes 
and Xc diverges at the phase boundary. Obtained CO 
pattern is schematically shown in the inset of Fig. 3(b). 
This is a ferri-electric structure and is similar to the one 
proposed in Ref.^^-' Within a chain along the x axis, the 
dipole moments align uniformly. Directions of the dipole 
moments in sublattices A and B are almost perpendicu- 
lar with each other. As a result, the electric polarization 
appears along the y axis. This ordered pattern is ener- 
getically favorable for the intcr-sitc Coulomb interaction 
Wij and the term which is proportional to QfQ^j in T-Lj. 

We also examine the dielectric and magnetic struc- 
tures in case II by the MC method. Temperature de- 
pendences of the charge susceptibilities Txc are pre- 
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Fig. 4. (Color online) Temperature dependences of the charge 
susceptibility Txc obtained by the MC simulation in the case 
a) . The inset shows the susceptibilities devided by L in the case 
c). The spin susceptibility Txg is also plotted. 

sented in Fig. 4 for several t/to. It is shown that with 
increasing L, Txc^s in t/to ~ 0.88 and 0.92 tend to 
be constants. On the other hand, the susceptibihty di- 
vided by the system size, Txc/L, in t/to = 1 is almost 
independent of L for T/to < 0-05- Therefore, the CO 
phase realizes in T/to < 0.05 in t/to = 1; f^nd oth- 
ers are the DM phase. The results in t/to — 0.88 and 
0.92, and the ones in t/to = 1 correspond to the case a) 
and the case c), respectively, although the corresponding 
values of t/to Si,iid the ordering temperatures are differ- 
ent from those in the MF calculations. With decreas- 
ing T in t/to = 0.92, Txc increases and suddenly goes 
down around T/to — 0.03 where the spin susceptibility 
Xs[= (TL)-iE,^.e'i--M(S,.S,)] at q = (1/3, 1/3), corre- 
spending to the 120° structure, develops. The results are 
consistent qualitatively with the MF calculation results. 

Finally, we discuss implications of the present results 
for the dielectric and magnetic properties in k-(BEDT- 
TTF)2Cu2(CN)3. We at first speculate that the transi- 
tion temperature for the 120° spin structure seen in Fig. 2 
(b) is interpreted to a temperature where the AFM spin 
correlation is developed in this compound. This is rea- 
sonable since this spin structure appears in a DM phase 
where spin frustration survives and the classical order- 
ing temperature is expected to be largely reduced. On 
the other hand, calculated spin orders in the CO phase 
and their transition temperatures are substantial, be- 
cause the anisotropic exchange interactions in the CO 
pattern shown in Fig. 3(b) release spin frustration. 

Based on the present study, we propose two 
possible scenarios which are relevant to k-(BEDT- 
TTF)2Cu2(CN)3. In the first scenario, the dipole mo- 
ments freeze in low temperatures. This corresponds to 
the sequential phase change of the case c) introduced 
previously; (Q:DM, S:Para)-> (Q:CO+DM, S:Para) ^ 
(Q:CO-|-DM, S:6-fold). A relaxor-like dielectric disper- 
sion observed experimentally requires a random freezing 
or a dipolar glass state below temperature where the di- 



electric constant takes a peak. This may be consistent 
with the experimental results suggesting inhomogeneous 
magnetic moment in low temperatures.'^^' We propose 
from the present calculations that the characteristic mag- 
netic temperature and the spin correlation functions in 
the CO phase are smaller than those in the DM phase, 
because the CO and magnetic phases are exclusive with 
each other as shown previously. One problem in this sce- 
nario is that the random/uniform polarization freezing 
makes the exchange interaction strongly anisotropy and 
releases spin frustration. This may be unfavorable for 
the experimental fact that neither a long-range mag- 
netic order nor a spin glass transition appears down to 
32mK.i"'") 

Another possible scenario is that the dipole moments 
do not freeze in low temperatures, and the dielectric 
anomaly experimentally observed is caused by charge 
fluctuation at vicinity of the phase boundary. This cor- 
responds to the sequential phase change of the case a); 
(Q:DM, S:Para)-> (Q:DM, S:120). Because of the reen- 
trant feature of the DM phase shown in Fig. 2(b), the 
system approaches to the phase boundary at certain tem- 
peratures, and the charge susceptibility is enhanced as 
shown in Fig. 3(a) and Fig. 4. With decreasing temper- 
ature furthermore, development of the spin correlation 
stabilizes the DM phase and suppresses the dielectric 
fluctuation. In this scenario, spin frustration is alive in 
the low temperature DM phase. Further examinations 
will provide a clue to unveil the microscopic pictures in 
the low temperature magnetic and dielectric structures. 
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